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FOREWORD 


The  Betasonde  was  originally  developed  by  Panametrics,  Inc.  under 
the  sponsorship  of  the  Division  of  Isotopes  Development  of  the  U.  S.  A.  E.  C. 
(Contract  Monitor:  Mr.  John  C.  Dempsey)  and  the  Atmospheric  Sciences 
Branch  of  O.  N.  R.  (Scientific  Officer:  Mr.  James  Hughes).  Several  launches 
of  these  early  versions  of  the  Betasonde  were  carried  out  at  While  Sands 
Missile  Range,  under  the  supervision  of  Mr.  Harold  N.  Ballard  of  the  Atmo¬ 
spheric  Sciences  Laboratory. 


The  advanced  model  Betasonde  work  reported  here  was  supported  by 
the  Army  Research  Office  and  was  monitored  by  Dr.  Arthur  V.  Dodd  of 
the  Geosciences  Division.  The  objective  of  the  work  was  to  develop  a  flight 
model  of  the  Betasonde  for  integration  into  an  Areas  rocket  payload  and 
application  in  a  middle  atmosphere  flight  program  al  While  Sands  Missile 
Range.  This  objective  has  been  accomplished  and  the  sonde  is  now  avail¬ 
able.  However,  the  WSMR  program  in  which  the  sonde  was  to  lie  used  was 
discontinued  during  Hie  fabrication  phase.  Should  funding  be  available  for 
a  flight  program  of  density  measurements  up  to  100  km,  the  Betasonde 
could  be  integrated  into  any  of  a  variety  of  rockets  and  used  for  the  purpose 
int  ended. 
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INTRODUCTION 


1.1  Betasonde  Development  Summary 

The  importance  of  air  density  measurements  has  long  been  recognized. 
Furthermore,  the  importance  is  increasing  due  to  the  necessity  of  observing 
and  forecasting  atmospheric  variables  for  re-entry,  SST  and  other  high  per¬ 
formance  vehicles.  Under  sponsorship  of  several  Government  Agencies, 
Panametrics,  Inc.  developed  an  Areas  borne  sonde  that  utilizes  beta  ray  for¬ 
ward  scattering  for  density  measurement  up  to  about  60  km  as  shown  below. 


Figure  1.  1  Cross  Section  of  Original  Betasonde  Showing  Density 
Measurement  Concept 


Essentially,  a  shield  is  interposed  between  an  annular  shaped  beta  source 
and  a  detector  in  such  a  manner  that  no  direct-path  beta  particles  can  be 
detected.  In  the  high  altitude  region  (>  20  km),  the  detector  counting  rate 
is  then  almost  linearly  dependent  on  the  atmospheric  density,  since  only 
air-scattered  particles  can  reach  the  detector. 

Under  AEC,  ONR,  WSMR  and  ARO  sponsorship,  there  have  been  three 
models  of  this  device  designed  and  fabricated: 

1)  Original  Betasonde 

2)  Low  Background  Betasonde  (LBB) 

3)  Betasonde  II. 


The  work  is  detailed  in  References  1.1-1.  3.  The  original  Betasonde 
used  standard  size  geiger  counters  (which  had  a  substantial  cosmic  ray- 
background  effect),  whereas  the  LBB  used  miniature  ones  (the  same  type 
used  in  discovering  the  Van  Allen  Belts)  to  reduce  this  effect.  Betasonde  II 
uses  a  semiconductor  detector,  which  has  a  very  low  cosmic  ray  background 
effect.  It  is  the  subject  of  the  effort  described  in  the  present  report. 

All  three  were  designed  for  flight  on  an  Areas  rocket,  and  when  inte¬ 
grated  into  it,  appear  essentially  as  above.  The  density  measurements 
are  made  as  the  sonde  descends  by  parachute;  hence,  it  is  recoverable  and 
reusable.  As  the  work  has  progressed,  we  have  pushed  the  upper  altitude 
limit  higher  with  each  model;  Betasonde  II  could  require  use  of  a  boosted 
Areas  with  an  apogee  in  the  90  -  100  km  region,  provided  a  sufficiently  in¬ 
tense  beta  source  is  used.  Results  obtained  with  the  original  version  are 
discussed  in  Ref.  1.  3. 

Modification  of  the  original  betasonde  to  use  the  miniature  geiger 
counters  was  funded  by  White  Sands  Missile  Range  (WSMR)  and  ONR.  It 
was  flown  very  successfully  to  an  altitude  of  about  65  km  at  WSMR  on 
April  23,  1976,  and  a  report  (PANA-AIR-1,  Ref.  1.2)  is  available.  Those 
atmospheric  density  measurements  are  included  in  Ref.  1.  1  (Figures  3.  1 
and  3.  2,  pg.  23  -  24). 

During  the  development  work  on  the  LBB  we  became  convinced  that 
the  technique  could  be  extended  to  the  90  -  100  km  region  by  use  of  a  semi¬ 
conductor  detector  and  a  more  intense  beta  source  («500  mCi).  Following 
submission  of  a  proposal  to  the  ARO  in  November,  1975,  wc  began  work 
on  a  planned  three-phase  Betasonde  II  effort: 

Phase  I  -  Research  and  Development 

Phase  IF  -  Design  and  Fabrication 

Phase  III  -  Flight  and  Data  Analysis 

The  ARO  Rinding  was  for  the  first  two  pha  •«,  with  WSMR  expected  lo 
fund  the  third  phase.  Near  the  end  of  the  second  phase,  however,  we  were 
informed  that  sucli  atmospheric  research  can  no  longer  be  funded  by  WSMR. 

Our  progress  througii  the  first  phase  is  detailed  in  the  annual  report, 
Ref.  1.  1.  That  work  is  briefly  summarized  in  the  present  work,  along 
with  details  of  the  design  and  fabrication  work  in  Phase  II.  That  effort 
has  now  been  completed,  and  we  could  undertake  Phase  III  work  at  any  time. 

That  work  could  be  described  as  follows:  the  instrument  must  be 
integrated  into  an  Areas  payload  for  launch  at  Wliite  Sands  Missile  Range 
(an  alternate  vehicle  or  range  could,  of  course,  be  used),  calibrated  at 


both  Panametries  and  in  a  large  (40'  -  60')  chamber  in  order  to  verify  the 
high  altitude  (>  70  km)  portion  of  the  calibration  curve.  The  proper  inten¬ 
sity  beta  source  must  be  procured  for  this  work.  The  Retasonde  II  should 
then  be  flown  to  the  80  -  90  km  region  at  least  twice,  and  the  analytical 
fitting  procedure  described  in  PANA-AIR-2  (Ref.  1.  I)  should  be  applied 
to  the  data.  This  will  form  the  basis  for  routine  direct  measurement  of 
atmospheric  density  up  to  at  least  the  80  -  90  km  region.  Reports  would 
be  prepared  and  the  results  published. 

We  would  like  to  point  out  that,  duo  to  the  development  of  lower  noise 
semiconductor  detectors,  the  possibility  exists  that  this  technique  could, 
in  fact,  be  extended  to  the  region  slightly  above  100  km.  Rut  we  believe 
that  such  a  further  development,  if  desirable,  should  take  place  only  after 
testing  Retasonde  II  and  determining  its  upper  altitude  limits. 

1.2  Detailed  Description  of  Retasonde  II  Work 

The  Scope-of-Work  for  the  subject  contract  is  as  follows: 

The  Cont  factor,  as  an  independent  Contractor  and  not  as  an  agent  of 
the  Cove  rnment ,  shall  provide  the  necessary  management,  facilities,  ser¬ 
vices  and  materials,  as  required,  to  accomplish  the  research  project  en¬ 
titled  "ATMOSPHERIC  DENSITY  MEASUREMENT  IN  THE  MIDDLE  ATMOS¬ 
PHERE’’  in  accordance  with  the  Contractor's  proposal  dated  November  1978 
and  revised  during  April  1976  which  is  incorporated  as  part  of  this  contract 
I > v  reference.  The  Contractor  shall  use  the  types  of  personnel  listed  in  Ex¬ 
hibit  A,  attached,  at  approximately  the  level  of  effort  stated.  The  research 
to  be  conducted  will  include  but  will  not  necessarily  he  limited  to  the  design 
and  fabrication  of  an  instrument  utilizing  beta  radiation  to  measure  atmos¬ 
pheric  density  t  o  a  n  c  lc\  a  t  i  on  ol  10)  kilometers.  The  1  i  r  s  t  pi  ,a  se  <  0  the  re  - 
search  involves  investig.'lt  ion  of  satisfactory  beta  radiation  sources  and  de¬ 
tection.  Space  simulation  chandlers  will  he  used  to  measure  and  i  alculate 
the  density  response  to  beta  radiation.  Rased  on  the  results,  a  densits 
sonde  w  1 11  he  designed  a  nd  lahricatrd  in  the  second  phase  ot  the  work. 

[bus,  t  ho  work  under  this  contract  was  divided  into  two  phases,  as  de¬ 
scribed  above,  and  both  have  been  sxicc.es sfu! lv  completed. 

Section  2  contains  a  description  of  the  method  of  density  measurement 
by  use  of  beta-ray  forward  scattering.  A  technique  is  presented  for  uti¬ 
lizing  an  entire  array  of  density  measurements  to  determine  the  optimum 
fit  to  an  analytic  altitude  dependent  profile.  Measurements  made  with 
the  Low  Background  Retasonde  are  used  to  demonstrate  the  technique. 


The  design  and  fabr icat inn  work  on  Betasonde  II  arc  detailed  in 
Section  3.  As  discussed  in  Sections  2  and  3,  calibration  results  for 
Betasonde  II  yield  similar  density  sensitivity  to  that  for  the  LBB,  pro¬ 
vided  the  same  radiation  source  is  used.  For  atmospheric  application, 
the  principal  difference  between  the  two  sondes  is  that  Betasonde  II  has 
very  low  cosmic  ray  background  sensitivity,  and  so  has  a  higher  ulti¬ 
mate  altitude  limit.  To  apply  Betasonde  II  up  to  the  80  km  region,  it 
will  be  necessary  to  obtain  the  proper  beta  source  having  strength  of 
about  500  mCi.  The  basic  specifications  for  Betasonde  II  are  given 
at  the  close  of  Section  3. 

Sectii  n  4  contains  a  short  summary  of  results  accomplished  in  this 
program  and  the  conclusions  regarding  potential  application  of  Betasonde  II. 

2.  ATMOSPHERIC  DENSITY  MEASUREMENT  METHOD  AND  RESULTS 

In  this  section  the  basic  concept  of  measuring  atmospheric  density 
by  use  of  beta-ray  forward  scattering  is  summarized,  and  the  accuracy 
that  can  be  expected  in  any  single  measurement  is  presented  based,  essen¬ 
tially,  on  the  statistics  of  the  counting  process.  In  any  actual  flight  of 
the  sonde,  there  would  be  a  long  series  of  such  measurements.  Thus,  re¬ 
sults  are  next  presented  for  a  method  that  uses  all  of  the  measurements 
to  find  the  best  fit,  in  a  least  squares  sense,  to  the  usual  exponential  de¬ 
pendence  of  density  on  altitude.  All  of  tire  2-csults  in  this  section  have 
been  given  previously  in  more  detail  (Ref.  1.  1,  Sections  I  and  4). 

2.  1  Description  of  Method 

The  measurement  of  atmospheric  density  by  beta-ray  forward 
scattering  relies  upon  the  fact  that,  at  low  densities,  the  number  of  par¬ 
ticles  scattered  by  a  gas  surrounding  a  beta  source  is  proportional  to 
the  density.  The  source  and  detector  must  be  arranged  in  a  geometry 
such  that  the  particles  can  only  be  detected  if  they  scatter  from  the;  gas, 
as  in  the  forward  scattering  configuration  shown  in  Figures  1.  1  and  2.  1. 

The  single- scattering  count  rate  C_s(Ej.)  cps  of  the  detector  for  the'  source' 
s  and  detector  energy  threshold  Ef  (the  energy  above  which  all  electrons 
are  counted)  is  found  by  integrating  over  all  scattering  angles  0  into  solid 
angles  d$2: 

Cs(K()  p(nNo/M)yy*  S  (E(  UMEtd.HE.  OklE  (2.1) 

il  E 

t 

Here  atmospheric  density,  s/cm' 

N  Avogadro's  number,  (>.  02  1x10^  ’  mol  ec  ulos/:  ;■  -  mol  ■■ ) 
o 

M  molecular  weight  of  atmosphere,  g-molc 

n  number  of  atmospheric  atoms/ mol  ccul  e 

dir  f(t))Et  EldSl  cm^/atom 


f 


FIE)  (l-(i2)/p4 

i(0)  :  (Z^rt)2/4)  Sin'^y/2  cm“/(»r-atom) 
p  -  v/c  -  ratio  of  beta  particle  velocity  to  that  of  light 
Z  atomic  number  of  atmosphere 
rQ  classical  electron  radius,  2.82x10“*  3cm 

Ss  lEj)  -  total  emission  rate  of  source  for  particles  of  energy 
greater  than  Ef,  jl/sec 

R  (E)  -  spectrum  for  beta  source  of  type  s.  normalised  to  unity 
from  Et  to  -x  ,  |1  /  k  e  V 
E  kinetic  energy,  keV 

The  cross  section  dct(E,  0)  for  scattering  of  an  electron  of  energy  E  (in  the 
energy  range  of  interest  here,  >50  keV)  through  an  angle  0  is  given  by  Evans 
(Ref,  2.  1).  As  shown  above,  it  is  separable  into  two  functions,  the  first  de¬ 
pending  onlv  on  E,  and  the  second  only  on  0 .  Thus,  as  shown,  the  count 

rate  is  proportional  to  p,  and 


C  ,(E1  =  K  (E  )()  cps 

St  Si 


(2.  2) 


Clearly,  K  depends  only  on  the  source,  detector  and  geometry  charac¬ 
teristics  as  given  by  the  right  hand  side  of  Eq  ,  (2.  1),  It  can  be  determined 
analytically  for  any  particular  beta  source  and  geometry.  This  has  been 
done  (Ref.  1.1,  Section  4)  with  the  result  that  for  those  sources  presently 
available,  Pm-147  is  the  best  choice.  It  has  a  suitably  long  half  life, 

2.  5  years,  and  a  maximum  energy  of  225  keV -which  is  sufficiently  high 
that  about  half  of  t  lie  emitted  beta  particles  have  energy  in  excess  of  the 

50  keV  threshold  E  (See  Ref.,  2.2  for  beta  spectra).  'I'hu.s,  for  the 
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Figure  2.  1  Reta-ray  Scattering  Configuration 


present  system,  we  believe  that  Pm-147  is  optimum.  It  will  provide  a 
means  of  establishing  the  utility  of  the  betasonde  for  the  upper  portion 
of  the  middle  atmosphere.  Results  obtained  by  flight  of  a  relatively  low 
intensity  Pm-147  source  on  the  Low  Background  Betasonde  (LBB)  can 
be  used  to  estimate  the  accuracy  that  would  be  obtained  with  higher  in¬ 
tensity  sources.  Before  doing  this,  it  is  necessary  to  discuss  the  experi¬ 
ment  al  calibration  procedure  by  which  K  can  be  determined  for  a  particular 
sonde  and  source. 

A  difference  between  the  present  sonde  and  the  LBB  is  that  a  signif¬ 
icant  wall  effect  is  present  during  calibration  in  the  Panamctrics  chamber 
(about  I'xV).  This  was  not  experienced  by  the  LBB  duo,  we  believe,  to 
tile  extremely  collimated  geometry  of  the  geiger  counters  themselves. 

Tims,  going  to  the  semiconductor  letectors  causes  a  significant  increase 
in  overall  efficiency,  hut  also  produces  a  wall  effect  in  the  small  chamber. 
This  does  not  affect  the  determination  of  the  constant  Ks,  since  it  is  the 
slope  of  the  count  rate  versus  density  curve.  It  does  mean,  however,  t  ha  t 
in  order  to  obtain  an  actual  calibration  curve  at  densities  less  than  that 
equivalent  to  about  70  km,  it  will  be  necessary  to  use  a  much  larger  cham¬ 
ber.  Such  chambers  are  available  at  NASA  Langley  Research  Center  (both 
’O'  ;ind  o0  diameter  spheres)  and  have  been  used  previously  (Ref.  J.  I)  for 
belasonde  calibration  work. 

Once  it  is  established  that  a  calibration  measurement  Ks  in  a  small 
chamber  leads  to  the  same  result  as  that  obtained  in  a  large  chamber  -  in 
which  the  wall  effect  is  much  reduced  -  it  is  not  then  necessary  to  recalib¬ 
rate  in  the  large  chamber  unless  the  source-detector  geometry  is  changed 
in  some  significant  way.  Rather,  a  value  of  Ks  can  be  carefully  measured 
in  tile  small  chamber  in  the  /,  <  t> 0  km  region,  where  the  wall  effect  is  only 
a  fraction  of  total  count  rate, even  in  the  small  chamber,  and  the  <  alibration 
curve  for  higher  altitudes  can  he  obtained  by  extrapolation.  This  is  possible 
because  single-scattering  is  by  far  the  dominant  effect  for  z  >40  km,  hence 
the  calibration  curve  is  extremely  linear  and  the  count  rate  versus  density 
is  riven  rigorously  by  Eq.  (2.2). 

Thus,  the  calibration  results  in  the  small  chamber  c.m  be  used  to 
estimate  the  accuracy  obtainable  at  high  altitudes  and  with  a  more  intense 
sou  r  e  e. 


The  source  strength  of  the  same  annular  Pm-147  beta  sourc  e  p;  c- 
viou:  !y  flown  on  the  Low1  Background  Betasonde  (LBB  I  was  measured 
on  M.i\  12,  1  ' ’7 7  to  he 


2.  7f>  x  1  0  '  |V  see 


U 


Ss(E,l 


7  4.’'  m( ’  i 

(lest  sour 


for  =  50  keV.  The  slope  of  the  count  rate  versus  density,  Eq.  (2.  2),  was 
measured  in  the  geometry  of  Fig.  1.  1  to  obtain  the  calibration  constant 

Kg(E()  =  0.93  x  i0<j  cps/(g/cm3)  (2.4) 

(test  source) 


for  the  source  intensity  (2.  3).  The  efficiency  of  Beta  sonde  II  was  found 
to  be  somewhat  better  than  that  of  the  L13B,  but  for  purposes  of  illustration, 
we  use  the  LBB  calibration  results. 


Now,  the.  fractional  statistical  uncertainty  for  any  count  rate  C  measured 
for  a  time  t(sec)  is  given  by  (Ref.  1.2) 


Ap/p  --  l/\/c  t 


(2.  5) 


By  use  of  (2.  2).  this  is 


Ap/p  =l//ks(Et)pt 


(2.  6) 


This  shows  explicitly  that  the  statistical  uncertainty  is  a  function  of  the 
threshold  E(  through  the  calibration  constant  KS(E^).  Of  course,  Ks  is  also 
proportional  to  the  source  strength  Ss(E^)  through  (2.  1),  We  believe'  that  a 
source  of  intensity  about  5  times  higher  than  (2.  3)  should  be  obtained  to  use 
<i  s  (lie  first  flight  source, 


S  ( E  )  =  1 . 38  x  1  01  ^  |l/s  ec:  •-  37  3  mCi, 
s  t 

(flight  source) 


for  which 


o  3 

Kg(E^)  r  4.h5  x  10'  cps/(g/ent  ). 

(flight  sour  t  e) 


(2.  7) 


(2.  8) 
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Now,  by  list-  of  the  standard  atmosphere  densities  (Ref.  2.  3),  Eqs.  (2.2), 
(2.  8)  and  (2.  6),  the  following  results  are  obtained: 


Table  2.  1 

Estimated  Fractional  Statistical  Uncertainty  for  Betasonde  II 
and  First  Flight  Source  (S^~  37  3  mCi) 


/ 

c 

Ap/p 

P  3 

km 

g/cnr 

cps 

t  =  20  sec 

t  -  5  sec 

t  =  1  sec 

60 

3.1x1  O'7 

1441 

.006 

.012 

.026 

70 

8. 8x1 O'  8 

409.2 

.01  1 

.022 

.049 

80 

2.0x1  O'8 

93.0 

.02  3 

.046 

.104 

0() 

3.2x1  0'9 

14,0 

.058 

.116 

.259 

1  00 

5. 0x1 O'l  0 

2.3 

.147 

.295 

.659 

The  uncertainties  listed  in  this  table  are,  of  course',  the  minimum  t  hat  will 
occur  for  the  indicated  count  rates.  The  actual  total  error  will  be  slightly 
larger  due  to  some  (as  yet  undetermined)  cosmic  ray  effect.  This  is  ex¬ 
pect  ed  to  hi'  extremely  small,  however;  much  less  than  in  the  I. HI’,  (Hef. 
1.2). 


A  typical  rocket  payload  spends  about  20  seconds  within  a  1  km 
region  near  apogee.  As  noted  above,  this  allows  a  very  accurate  mea¬ 
surement  of  the  density  in  that  particular  region.  But  as  the  payload 
descends,  it  initially  gains  speed  until  it  is  slowed  by  entry  or,  if  the 
payload  is  to  be  recovered  as  for  an  Areas,  by  a  parachute.  There  is, 
then,  a  region  in  which  the  speed  is  high  between  apogee  and  entry. 

In  this  region  the  individual  measurements,  if  assigned  to  altitude  in¬ 
crements  similar  to  those  at  higher  and  lower  altitudes,  are  less  accurate. 
Thus,  a  procedure  has  been  developed  that  utilizes  all  of  these  measure¬ 
ments,  with  their  associated  accuracies,  to  define  a  density  profile. 

This  profile  is,  effectively,  "tied"  at  high  and  low  altitudes  to  some  very 
accurate  measurements,  with  the  mid-altitude  portion  required  to  fit 
the  less  accurate  data  in  that  region  and  to  fit  smoothly  with  the  profile 
for  the  other  regions. 

In  the  following  sub-section,  this  procedure  is  summarized,  and 
it  is  applied  to  results  for  a  1,1313  flight  to  demonstrate  its  applicability. 
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2.  2  Opt iirmm  Analytical  Fit  to  Altitude- Density  Profile 
2.  2.  1  General  Objective 


As  is  well  known,  density  ()(/.)  is  approximately  an  exponent  jal  I  y 
decreasing  function  ot  altitude 


f>(z)  -  Ps1' 


-z/H 


where  H  is  the  scale  height  and  ps  the  value  of  p  at  ■/,  0.  In  fact,  stub  an 
approximation  is  often  made  use  of  (Ref.  2.  -1,  p  o7)  to  provide  analytical 
rep  resenlations  of  models,  for  use  in  applications  win- re  such  an  approach 
is  more  convenient  than  is  a  long  tabulation,  for  example  whi  rr  a  i  ampu¬ 
tation  is  to  be  made  that  requires  a  continuous  representation  versus  alti¬ 
tude.  The  result  is  generally  provided  in  the  form  of  the  coefficients  in  a 
power  series  expansion  lor  Pnp(y.).  he  believe  that  a  similar  approach 
could  also  be  useful  lor  representing  the  data  from  a  single  determination 
ot  a  (>(■/.)  profile,  in  particular,  from  a  rocket  launch  such  as  that  ot  Pu  la- 
sonde  II,  the  instrument  developed  under  the  present  contract  for  air 
density  measurement  up  to  at  least  80  km  by  beta  ray  torward  scatterin'.  . 

1  'he  1  ,ow  -  Rac  kg  round  Uetasunde  (I.HH)  is  an  earlier  version  Areas 
borne  instrument  designed,  prineipal  1  y ,  for  measurement  ot  i >{•/.)  between 
about  28  and  to  km.  The  direct  output  ot  tiie  data  reduction  procedure  is 
a  tabulation  of  individual  densities  p ■(•/.)  and  associated  statistical  uncer¬ 
tainties,  Apj(z),  versus  ■/.-  When  making  use  of  data  of  this  type  it  is 
normally  not  possible-  simply  U>  interpolate  directly  between  values  in  tin- 
tabulation.  The  reason  is  that  tin-  points,  when  plotted,  will  generally 
show  an  amount  of  scatter  consistent  with  the  uncertainties.  Tims,  in 
order  to  use  such  data  it  is  necessary  to  draw  a  best-fit  curve  of  some 
type  nil  the  plot  (usually  simply  by  eye).  Values  are  then  read  oil  the 
curve  and  assumed  to  have  an  uncertainty  consistent  with  that  of  tin-  sta¬ 
tistical  Ap'(z),  and  any  additional  systematic  errors,  in  that  altitude  re¬ 
gion.  Clearly,  then,  determination  of  an  optimum  analytical  lit  lor  a  tab¬ 
ulation  would  serve  two  purposes.  It  would  provide  a  ''formal'1  method  ot 
lilting  the  data,  and  the  resulting  analytical  fit  would  provide  a  convenient 
representation  lor  computational  purposes. 

We  will  use  the  results  of  a  launch  "f  the  1,1111  at  White  Sands  Mi  - 
sile  Range  at  I  7  ;0  0  :0  0  l!T  on  April  2 t ,  l'.)?b  (Kef.  1.2)  to  show  the  i  y«- 

analytical  fit  that  ran  he  obtained.  First  the  general  matliematii.il  ap¬ 
proach  is  presented,  and  this  is  followed  bv  tile  densitv  ta  In:  hit  ion  ,o  :  -he 
ipi-lii.it  ion  ot  the  approach  to  the  tabulation  in  determine  ihe  ‘it. 


2.  2.  2  Mathematical  Approach 


The  deviation  oi'  (2.  9)  from  on  exact  exponential  <.  an  he  taken  into 
account  by  allowing  the  scale  height  H  to  he  a  lunction  of  altitude.  Over 
the  altitude  region  up  to  about  100  km,  11  is  not  stroii-;])  dependent  m  z 
and  should  be  representable  as  a  simple  iunction  (this  is  essentially  what 
is  done  in  the  standard  approaeh  to  model  fitting).  Ii\  doing  this  the  re  - 
suit  H(z)  is  a  physically  meaningful  parameter,  although  lor  a  given  anul- 
vtieal  fit  the  dependence  on  z  can  he  expected  to  ag  rue  wit!)  model  re  cuts 
in  onl\  an  approximate  manner.  For  purposes  ot  determining  a  :i!  ser  a 
particular  altitude  region  it  is  then  convenient,  and  appropriate,  t  >  ■  use 


P 


o 


z(1)/H(z) 


12.  1  0) 


when'  z  is  the  (specified)  minimum  altitude  of  interest.  i<(,  is  the  value 
iif  p  at  z  zu  and  ll(z.)  is  an  "average"  value  of  !l(z),  as  determined  by  the 
fitting  procedure.  Thus,  once  pG  and  lT(z.)  are  found  from  the  procedure, 
the  analytic  al  r e.sult(2.  ]  0)r (-presents  the  entire-  array  of  measured  density 
profih-  data,  which  are  in  the  form  and  statistical  uncertainty  for  all 
N  altitudes  z.-.  The  procedure  is  used  to  minimize  tlu  sum  of  the  mean 
squared  fractional  deviations  Dj  of  all  p  values,  calculated  from  (2,  10)  from 
the  measured  results  pp 

l>.  (P  -pjJ/p  ^  1  -P-JP  1)1 


In  order  to  determine  how  well  the  lit.  has  been  made,  once  the  pro¬ 
cedure  is  complete,  it  is  useful  to  compare  two  quantities.  The  first  is 
the  rms  uncertainty  Am  of  the  N  values  of  I  raelioual  unci  rtainly  Aj  in  the 
i  ii. i i \  i d u a  1  measurement, 


A 


N 


/.\ 


(2.  1  .•> 


W'll  l*  IV 


A.  Ap./p. 

i  \  '  i 


</.  1 


.  i  iw  1  ,\|)  j  is  tlu*  statistical  unci-  rl  a  ini  y  l^ci.  !.*-»  i>  I  '">  )• 
i  s 


A 


> 

c 


Z(Pi/p  •■'■0  -  ')7>- 


Tile  <  rum:  i  p  la  lit  it  y 


1  0 


!  2.  ill 


which  is;  the  rms  deviation  of  (lit-  p(zj)  v .  t  J  i  ( t  ■  s  .  e.iic.  ini.  <|  by  usi-  >,i  [2.  10), 

1  ivm  the  mo  a  s  u  rad  values  p^.  Ii  <ili  moa  s  u  n.  n  leni ■./<  re  mark-  a)  only  <■  no 
vaJuo  of  p,  then  it  would  bo  exported  that  A  ~A  ..  ' i  1  •  i is  ijn  .mn-,  f.,r 

given  actual  p,  the  random  variations  in  the  measured  p-  for  tl.,-  1'ot.tsonde 
are  duo  to  the  statistical  nature  of  the  nuclear  decay  pro<  ess  that  leads  t  > 
beta  emission.  The  values  of  p^  are  distributed  in  a  1  'oisson- like  manner, 
for  which  the  fractional  standard  deviation  of  the  mean  of  a  series  o!  N 
readings  would  he  given  by  (2.  12).  The  Pearson  "chi  -  qua  red"  lest  ear. 
readily  be  applied  in  such  circumstances  to  determine  how  well  the  mea¬ 
surements  fit  the  dist  ri  bution  as  sumed.  Here,  howe\er,  wi;  are  dealing 
with  the  comparison  of  u  best  analytical  lit  (rather  than  the  mean  o,  ...  series 
oi  measu  remolds)  to  the  individual  measurements,  each  <u  which  is  expected 
to  l)e  slightly  ditto  rent  because  /  is  changing.  Additionally,  sc, me  atmos¬ 
pheric  variability  will  exist  that  cannot  he  fitted  by  a  simple  Junction  such 
a  s  (  2. 1  0  )  without  an  unrealistically  large  amount  of  variation  in  H(z),  and 
this  will  in  itselt  cause  A  .  to  exceed  Ani.  Clearly,  in  tins  case  we  expect 
Ac  T-A.,,,  and  v.  c  can  conclude  that  the  ti!  is  sat  i  shie  t  o  r  y  i;  the  diifen  nee  is 
acceptable  small. 

Thus,  the  procedure  here  is  to  tinci  tne  simplest  possible  represen¬ 
tation  of  H(z)  that  will  produce  a  value,  of  A  that  is  a ecept.ioU  close  to  ,y 
A  lit  in  which  Al«Am  would  not  represent  a  true  best  t:  t ;  1 1 1 1 . ;  c  .add  he  ob¬ 
tained  by  using  a  sufficiently  com  plica  I  od  fum  lion  for  !!(/.)  urn1  would  c  o  r  - 
respond  to  following  the  random  deviations  in  the  me  a  u  r. -meats,  r.tth<-  r 
than  the  true'  a\  era  go  variation  of  p(z).  J  lie  question  o‘.  ,  >.,m  tlv  how  closely 
.\.  must  approach  Am  for  the  analytical  fit  to  repreceni  1 1 , ■  experimental 
tabulation  acceptably  is  not  considered  further  here.  Rather,  we  showed 
fl'ef.  1.2)  that  fits  can  be  obtained  in  which  the  difference  is  small  and 
reduces  in  the  manner  expected  as  the  number  of  A  increases. 

j 

For  the  fit  used  here,  a  set  of  coefficient  s  A.  was  determined  in 
which  we  represent  1I(/.)  in  an  inverse  power  series,  so  that 


f  Ml/p,)  A, 
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U  -  M/H(z.) 
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ind  j  is  the  number 
a\l  A  . 


of  coefficients  A. 


om« 

U.  U>): 
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are  known,  both  /)  and  H(zr)  can  b< 
o 


ciiosen  to  be  used, 
found  from 


Clearly. 

( 2 .  J  h)  and 


(2.  17) 


ll(z) 


/J">  ;  . 

7  >;.A "  -  *-'J " 


(2.  IK) 


Hence,  a  fit  must  he  made  to  determine  the  A.  for  use  in  (?,.  17)  and 
(2.  ]  .  ),  which  then  defines  p  ( z )  from  (2.  10)  ' 

i  Flight  Hat  a  and  Density  Tabulation 

Details  of  the  Dow  Background  JRtnsondc  (RUB)  flight  on  April  2  -1, 

I  D'u,  are  contained  in  Ref.  1.2  and  presented  here  only  to  the  e;-(--nt 
:'<•<.  t  ssary  for  application  of  the  above-  procedure. 

1  he  Areas  was  lnunclud  to  ati  apyer  altitude  near  n:  km,  with  the 
parachute  home  l.llil  bring  ejected  just  helore  apogee.  The  pa  r;<«  hole  ha., 
little  client  above  SO  km.  The  descent  velocity  changes  I  mm  about  km. 
sec  to  .02  km/sec  between  60  and  30  km.  However,  the  sense  does  pend 
tte.rly  21'  .seconds  within  1  km  of  apogee ,  thus  providing  much  data  is.  this 
region  before  beginning  its  rapid  descent.  This  is  at,  important  !..<  t, 
because  it  allows  the  density  profile  to  lie  "tied  down"  at  high  altitude. 
Clearly,  as  the  sonde  descent  rate  slows,  it  obtains  prog  re  s  si  v »  )  .  more 
data  in  each  altitude  region.  Hence,  the  entire  densitv  profile  for  the  low 
altitude  region  is  very  accurately  defined.  The  pmecdurt  deseriis  l  above 
is  designed  such  that  it  provides  a  single  a  t  • .  •.  i  y  t  i  t  a  1  til  that  automa  :  .e.i  11  y 
lakes  all  of  these  factors  into  account. 

The  BetasoJide  pulses,  each  of  whirl,  (except  :  r  iistim  rav 
g  round)  corresponds  to  detection  of  an  a  i  r  -  s,  a !  te  n  d  beta  part'  ,  r, 

te  I  i'iiii'I  e  red  amt  recorded  on  magnetic  lat'e  n.  reai-tiioi  .  1  he  .r  ia  . 

ore  tlnii  summed  into  approximalelv  eipia  i  i  i  i , .  *  ■  int  rv  I !  deoenDie 
ip. >n  telmiietrv  noise  and  time  ic'd  tor  transmitting  I  >  •  to  p>  ■  fa  ;  an  ■  ■  •  r  di- 
agnostic  purposes.  In  Table  2  these  sums  an  1..  be  1 1 1  i  t  "I  ,S  .  ■ :  >  ■  ,  h 
.  ( 1 1 1 1  tide  z  (1<  m)  I  ■  ■  r  t  In  lime  interval  '1  M(SK»\S).  1  he  e  ;i  n-  the  I...  :  .  !  : 

data  '  roil  i  who  h  tlie  densitv  m<  a  s  u  r,  meat  at  e.  .  b  .  ■  1 1  i !  u  I  e  is  ■  i  >  I  >  -  -  .  ■  •  .  ■  | , 


1  .1 


Table  2.2  Low  Background  Beta  sonde  Count  Bata  for  April  21,  l')7(, 


Z  (  KM  > 

TM< SECS) 

CTS 

64.0 

5.5 

80. 

64.3 

5.2 

89. 

64.2 

5.4 

85. 

63.5 

5.2 

83. 

62.9 

5.1 

102. 

62.1 

5.1 

99. 

61  .0 

5.2 

1  30. 

60.0 

5.2 

1  33. 

59. 1 

5.1 

156. 

57.9 

5. 1 

1  55. 

55.9 

5.  1 

219. 

54.9 

5.  1 

240. 

53.9 

5.  1 

259. 

53.0 

S.  1 

279  . 

52.0 

5.0 

355. 

51.2 

5.2 

376. 

50.5 

5.0 

376. 

50.0 

4.  ! 

311. 

49 . 5 

4.0 

361  . 

48.6 

5.0 

543. 

43 . 1 

5.0 

502. 

4  7.7 

5. 1 

565. 

47.2 

4.0 

467. 

4  6.8 

4.0 

48  5. 

4  6.4 

4.  1 

495. 

45.8 

2.0 

30  7. 

45. 3 

4.0 

650. 

45.0 

4 . 0 

626. 

44.2 

4.0 

683. 

43.9 

4.0 

748  . 

43.4 

5.0 

1047. 

43.  1 

5.  1 

989. 

42.8 

5.0 

1093. 

42.5 

5.0 

1  1  24. 

42.2 

5.  1 

1197. 

41  .9 

4.0 

942. 

41.7 

4.0 

1024. 

41.1 

5.0 

1  408  . 

40.8 

5.0 

1  443. 

40 . 5 

5.0 

1510. 

40 . 2 

5.  1 

1  5  38  . 

^0 . 0 

5.0 

1610. 

39.7 

5.0 

1  652. 

39.5 

5.0 

1  742. 

39. 3 

5.0 

1803. 

3^.0 

4.0 

1507. 

38 . 5 

4.  1 

1  574. 

38 . 3 

5.0 

200  7. 

33 . 1 

5.0 

208  6. 

37.9 

5.0 

2228  . 

3  7.6 

5.0 

2241  . 

37.4 

5.0 

2335. 

37.2 

5.0 

238  5. 

37.0 

5.0 

2400. 

36.8 

5.0 

2545. 

36.2 

5.0 

2808  . 

36.  ! 

5.0 

2964. 

35.9 

5.0 

2975. 

35.7 

5.0 

2890 

35.5 

5.0 

3102 

35. 4 

5.0 

3098 

35.2 

5.0 

31  66 

35. 1 

5.0 

3344 

34.9 

4.0 

2655 

34.5 

5.1 

3547 

34 . 4 

5.0 

3628 

34.  3 

5.0 

3530 

34.  1 

5.0 

3444 

34.0 

5.0 

3  6  6  6 

33.9 

5.0 

3809 

33.  7 

5.0 

38  29 

33*6 

5.0 

40  33 

33.5 

4.0 

3283. 

3  3.1 

5.0 

4097 

33.0 

5.  1 

4227 

32.9 

5.0 

4359 

32.8 

5.0 

43  6  3 

32.6 

5.0 

4580 

32.4 

5.0 

4735 

32.2 

5.0 

4755. 

32. 1 

5.0 

4895. 

32.0 

4.0 

3944. 

31  .  6 

5.0 

5055. 

31.5 

5.0 

5  353. 

31.3 

5.  1 

5  316. 

31  .  2 

5.0 

5271 

31  .0 

5.0 

5  406. 

30.8 

5.0 

5701 

30.7 

5.0 

5707. 

30.5 

5.0 

5901 

30 . 4 

4.0 

4  7  35. 

30.0 

5.0 

6  3  38 

29 . 9 

5.0 

6  418. 

29 . 7 

5.0 

6  3  32. 

29 . 6 

5.0 

6350. 

29.5 

5.0 

6  4  62. 

29 . 3 

5.0 

6  7  6  4 

29.1 

5.0 

6708 

29.0 

5.0 

6  7  5  7. 

28.9 

4.0 

5757. 

28 . 6 

5.0 

68  95. 

28 . 5 

5.  1 

7  38  3  • 

23 . 3 

5.0 

7  335. 

28.2 

5.0 

7  4  5  4. 

23 . 1 

5.0 

7  724. 

27.9 

5 . 0 

7  764. 

2  7 .  7 

5.0 

790  1  . 

2  7  .  t, 

5.0 

7  7  72  . 

1 . 4 

4 . 0 

6  5  02. 

2  7.2 

5.0 

8  4  ri  . 

2  7  .  1 

4.0 

6  3  4  1  . 

2  7.0 

■•■.0 

8  401  . 

36.9 

5.0 

8  2  34  . 

2  6.7 

5.0 

8  7  1  9  . 

2  6 .  6 

5.0 

8  6  75. 

2  6 . 5 

5.  1 

3  509. 

2  6.3 

5.0 

8914. 

2  6.2 

5.0 

8892. 

1  ' 


FI  i  L',ht 


it  is  soon  that  the  counts  gcm-rallj  im  rc;i:a-  .is  altitude  tl  ec  re.i  ses , 
because  tin'  number  of  beta  particles  scattered  increases  as  the  density 
increases  in  this  density  region.  In  addition  to  the  density  itself,  there 
are  several  other  factors  that  must  be  taken  into  ao  mint  in  determining 
density  from  this  count  rate  (see  Ref.  1.2  for  details  of  procedure): 

1)  circuit  dend-time,  caused  by  dead  time  in  i:eip  r  counters  after 
counting  a  pulse  (~25/.tsce)  and  a  (lead-time  circuit  introduced 
purposely  to  keep  the  count  rate  into  the  telemetry  transmitter 
from  exceeding  its  capabilities, 

2)  cosmic  ray  elfeet,  a  small,  but  non-n.gli  :ible,  altitude-depen¬ 
dent  e  Hi  cl, 

and  3)  time  elapsed  between  calibration  and  flight. 

K.iih  of  these  tactwrs  luts  ueett  taken  into  account  as  neeessarv  in  obtain¬ 
ing  the  calibration  data  of  Table  2.  3  and  the  i  light  dcnsii\  tabulation  in 
1  aide  2.  -f.  In  the  latter  table  the  count  data  1  rom  '1  able  2.2  arc  repeated 
as  the  lirsl  three  columns,  with  the  observed  counts  (O  I',. St  '.'l  S)  in  the  third 
This  leads  to  an  observed  count  rate  (not  sltowit)  which  is  corrected  mr 
dead  lime,  In  use  ot  the  circuit  de.nl  time  ot  1  .•i!'|j:;i'e,  to  obtain  the  t  orrec 
ted  count  rale  (COR.CTK.),  1  he  cosmic  r.iy  count  rate  (COSCl  It.)  is  then 
subtracted  fu  obtain  the  density-dependent  count  rate  (!)  KNC  1  R) . 

The  last  of  the  factors  that  must  be  taken  into  account  is  impor¬ 
tant  because  of  decay  in  intensity  of  the  radioact  ive  source,  in  this 
case  Pni-147  whose  mean  life  is  1  i 7 hi  days,  lass.  atially,  a  calibration 
curve,  that  would  bo  correct  on  the  day  of  the  flight  is  derived  from  that 
recorded  earlier.  In  the  present  instance,  the  tune  lapse  was  only  J  I, 
lays,  and  the  calibration  data  art^  shown  in  Table  2.  1  here  arc  22 

points  at  which  the  density  p(g/nT  )  and  the  count  rate,  corrected  for 
t  he  dead  time  factor,  was  as  shown  in  C  i  S/S  1:1  ( At  low  values  of  p, 
the  count  rate  increases  approximately  linearly,  but  as  p  increases 
the  count  rate  increases  less  rapidly.  At  2r  km  p  <10  g/m  and  at 
30  km  p  sr  20  g/m  .  Data  below  about  30  k.u  v,  •  •  r »  not  an  Tycd  in 
i\ef.  i.2  because  of  the  uon-Iinear  nature  of  the  eali.bi  at  Ion  curve  in 
lh.it  region.  Here,  however,  all  of  tin-  data  in  Table  2.2  are  ,ina!\ /, ed. 

This  was  done  by  making,  for  each  flight  count  rate,  a  quadratic  lit 
to  three  appropriate  calibration  count  rate  points  in  1'a'd  ■  2.  '.  II. e 
measured  PKNSITY  values  of  Table  2.  -I  were  tuns  determined  from 
dens  it  V  -  dependent  count  rate  DKNCTK,  and  the  analv'ieal  111  to  the  gen.  .ty 
count  rate  results  in  Table  2.  3. 

The  absolute  uncertainly,  UNC  (g/m  )  is  obtain,  d  from  :h<  density, 
the  counts,  and  the  calibration  data.  I  lie  fractional  uncertainty  2S  i:. 
given  by  the  ratio  IJNC/ DM  NS  IT  Y  . 
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Table  2.  3  Low  Background  Betasonde  Density  Calib¬ 
ration  Data  Recorded  April  7,  1976. 


PT  RHO<  G/M*+  3>  CTS/SEC 


01 

00 

0000 

02 

02 

0174 

03 

04 

0352 

04 

06 

0524 

05 

08 

0  700 

06 

10 

08  74 

07 

12 

1045 

08 

1  4 

1  220 

09 

1  6 

1  400 

10 

18 

1575 

1  1 

20 

1  750 

12 

22 

1  900 

1  3 

24 

2020 

1  4 

26 

21  30 

1  5 

28 

2245 

16 

30 

2345 

1  7 

32 

2440 

18 

34 

2530 

1  9 

36 

2620 

20 

38 

2705 

21 

40 

2790 

22 

42 

2880 

Table  2.  4  Low  Background  Beta;. en.de 
Tabulation  for  April  23,  1476  Blight 


Measu  red 


Density 


M  EAN-L I FE  *  DAYS- 1  378. 

TIME  FROM  CAL* DAYS- 16. 

CIK.  DEAD  TIME*SEC=. 000185 

UNITS:ALT=KM; TIME-INT=SECJ ALL  CTR  =  CPSi DEN* UNC=G/M** 3 


ALT  TIME-INT 

0BSCTS 

C0RCTK 

COSCTR 

DENCTR 

DENSITY 

UNC  FKAC-UNC 

64.0 

5.5 

80. 

1  4. 58 

0. 30 

1  4. 45 

0.168 

0.019 

0.112 

64.3 

5.2 

89. 

17.17 

0. 30 

1  7.07 

0.193 

0.021 

0.  106 

64.2 

5.4 

85. 

15.79 

0. 30 

15.67 

0.182 

0.020 

0.108 

63.5 

5.2 

83. 

16.01 

0.  30 

1  5.89 

0.185 

0.020 

0.110 

62.9 

5.  1 

102. 

20.0  7 

0.  30 

20.01 

0.232 

0.023 

0.099 

62.1 

5.  1 

99  . 

19.48 

0.30 

19. 41 

0.225 

0.023 

0.100 

61.0 

5.2 

1  30. 

25.  12 

0. 30 

25.11 

0.291 

0.026 

0.068 

60.0 

5.2 

1  33. 

25.70 

0. 30 

25. 70 

0.298 

0.026 

0.087 

59. 1 

5.1 

1  56. 

30. 76 

0. 30 

30.8  2 

0.  357 

0.029 

0.080 

57.9 

5.1 

1  55. 

30. 5  6 

0.30 

30.62 

0.355 

0.028 

0.030 

5  5.9 

5.  1 

219. 

43.29 

0. 30 

43.49 

0.504 

0.0  34 

0.067 

54.9 

5. 1 

240. 

47.47 

0. 30 

47.72 

0.553 

0.036 

0.064 

53.9 

5.  1 

259. 

51.27 

0.30 

51.56 

0.597 

0.037 

0.062 

53.0 

5.  1 

279. 

55.27 

0.30 

5  5. 61 

0.644 

0.0  38 

0.0  60 

52.0 

5.0 

355. 

71.94 

0. 30 

72. 48 

0.839 

0.044 

0.053 

51.2 

5.2 

376. 

73.29 

0. 30 

7  3 . 8  4 

0.R54 

0.044 

0.051 

50-5 

5.0 

376. 

76.26 

0. 30 

7  6.85 

0.889 

0 .046 

0.051 

1 


Tabk- 


4  (Coni . ) 


ALT 

T IME-INT 

OBSCTS 

COKCTK 

COSCTR 

denctk 

DEM  SI TY 

UNC  h  KAC-UNC 

50.0 

4.1 

31  1  . 

76.93 

0.30 

77.53 

0.897 

0.051 

0.056 

49.5 

4.0 

361  . 

91  .  78 

0.31 

92.54 

1  .069 

0.056 

0.052 

48.6 

5.0 

543. 

1  10.83 

0.32 

111.80 

1  .290 

0.055 

0.043 

48 .  I 

5.0 

502. 

102. 30 

0.32 

103. 1 7 

1.191 

0.053 

0.044 

47.7 

5.  1 

565. 

1  1  3.  10 

0.33 

11 4.09 

1  .3!  6 

0.055 

0.042 

47.2 

4.0 

467. 

119.33 

0.33 

1  20 . 38 

1  .  388 

0.064 

0 .  C  -16 

46.8 

4.0 

48  5. 

124.03 

0. 34 

125.14 

1  •  443 

0.065 

0.045 

46.4 

4.  1 

49  5. 

123.49 

0.34 

1 24. 58 

1 . 437 

0.064 

0.045 

45 . 8 

2.0 

30  7. 

1 57.99 

0.35 

159.48 

1  .835 

0.  1  04 

0.056 

45.3 

4.0 

650. 

1 67.54 

0.36 

169.13 

1.945 

0.075 

0.0  39 

45.0 

4.0 

626. 

161.17 

0.36 

1 62. 68 

1  .871 

0.074 

0.040 

44.2 

4.0 

68  3. 

1 76. 32 

0.37 

1 78.00 

2.044 

0.075 

0.037 

43.9 

4.0 

748. 

1  93.  70 

0.37 

195.58 

2.239 

0.079 

0.035 

4  3.4 

5.0 

1  047. 

217.84 

0. 33 

220.00 

2.510 

0.076 

0.0  30 

4  3.1 

5  •  1 

989. 

201.14 

0. 38 

203.10 

2.  322 

0.072 

0.031 

42.8 

5.0 

1093. 

227.81 

0.39 

230 .08 

2.622 

0.0  78 

0 . 0  30 

42. 5 

5.0 

1  1  24. 

234.55 

0. 39 

236.90 

2.699 

0.079 

0.029 

42.2 

5 . 1 

1  197. 

245.36 

0. 40 

247.82 

2.821 

0.080 

0.028 

41.9 

4.0 

9  42. 

246.23 

0. 40 

2  48 . 70 

2.831 

0.09! 

0.0  32 

41.7 

4.0 

1024. 

268 .73 

0.40 

271 . 46 

3.086 

0.095 

0.031 

41.1 

5.0 

1  408. 

297.08 

0.41 

300.1  3 

3.410 

0.091 

0.027 

40.8 

5.0 

1  443. 

304.88 

0.41 

308.02 

3.49  9 

0.092 

0.026 

40.5 

5.0 

1510. 

319.87 

0.42 

323.  1  9 

3.671 

0.095 

0.026 

40 . 2 

5.1 

1  5  38  . 

31  9 . 39 

0.42 

322. 69 

3.666 

0.094 

0.026 

40 . 0 

5*0 

1  610. 

342.40 

0.42 

345.97 

3.931 

0.098 

0.025 

39.7 

5.0 

1  652. 

351.91 

0 . 43 

355.59 

4.042 

0.103 

0.025 

39 . 5 

5.0 

1  742. 

372.40 

0.43 

376. 32 

4.285 

0.106 

0.025 

39.3 

5.0 

1803. 

38  6. 38 

0.43 

390.45 

4.451 

0.108 

0.024 

39.0 

4.0 

1  507. 

40  4.98 

0.43 

409.27 

4.671 

0.123 

0.026 

38.5 

4.  1 

1  574. 

41 3.25 

0.44 

417.63 

4.768 

0.123 

0.026 

38 . 3 

5.0 

2007. 

433.60 

0.44 

438.21 

5.008 

0.114 

0.023 

38 .  I 

5.0 

2086. 

452.09 

0.45 

456.92 

5.225 

0.116 

0.022 

37.9 

5.0 

2228. 

485. 63 

0.45 

490.85 

5.618 

0.120 

0.021 

37.6 

5.0 

2241  . 

488 . 72 

0.45 

493.97 

5.654 

0.120 

0.021 

3  7.4 

5.0 

2335. 

511.16 

0.45 

516.67 

5.916 

0.123 

0.021 

37.2 

5.0 

2  385. 

523. 1 7 

0.46 

528.81 

6.054 

0.122 

0.020 

37.0 

5.0 

2400. 

526.78 

0.46 

532.47 

6.096 

0.123 

0.020 

36.8 

5.0 

2545. 

561.91 

0 . 46 

568.01 

6.498 

0.128 

0.020 

36.2 

5.0 

2808. 

626.71 

0.47 

633.56 

7.242 

0. 1  36 

0.019 

36.1 

5.0 

2964. 

665.82 

0.47 

673.1 2 

7.693 

0.14! 

0.018 

35 . 9 

5.0 

2975. 

6  68. 60 

0.47 

675.93 

7.725 

0.141 

0.018 

35. 7 

5.0 

2890. 

647.21 

0  •  48 

654.28 

7. 478 

0. 1  39 

0.019 

35.5 

5.0 

3102. 

700 .8  4 

0  •  48 

708 .54 

8.097 

0.145 

0.018 

35 . 4 

5.0 

3098. 

699.82 

0. 48 

707.51 

8.086 

0.1  45 

0.018 

35 . 2 

5.0 

31  66. 

71 7.22 

0. 48 

725. 1 1 

8.286 

0.147 

0.018 

35.1 

5.0 

3344. 

763.23 

0. 48 

77} . 66 

8.819 

0.153 

0.017 

34.9 

4.0 

2655. 

756*  66 

0.49 

765.01 

8.743 

0. 1  70 

0.019 

34.5 

5.  1 

3547. 

798.19 

0.  49 

807.02 

9.226 

0.  1  56 

0.017 

34.  4 

5.0 

3623. 

8  38 . 1  0 

0.49 

8  4  7. 39 

9. 692 

0.163 

0.017 

34.3 

5.0 

35  30. 

8  12.06 

0.49 

821.05 

9. 388 

0.159 

0.01  7 

34.  1 

5.0 

3444. 

789 . 39 

0.50 

798 . 1  1 

9.123 

0.156 

0.017 

1  <> 


1  a!>lc  2.  4  (Coni.) 


ALT 

TIME-1NT 

OBSCTS 

COKCTK 

COSCTR 

DENCTk 

DEMSI TY 

UNC  FRAC-UNC 

34.0 

5.0 

3666. 

848.26 

0.50 

857.66 

9.81  1 

0.164 

0.017 

33.9 

5.0 

3809. 

886.78 

0.50 

896.63 

10.267 

0.171 

0.017 

33.7 

5.0 

3829  . 

892.20 

0.50 

902. 1 1 

10. 332 

0.172 

0.017 

33.6 

5.0 

40  38  . 

9  49 . 45 

0.50 

960.04 

11.012 

0.177 

0.016 

33.5 

4.0 

328  3. 

967.63 

C.50 

978.47 

11.227 

0.199 

0.018 

33.  1 

5.0 

4097. 

965.81 

0.51 

976.57 

11.205 

0.178 

0.016 

33.0 

5.1 

4227. 

973 .92 

0.51 

989.84 

1 1 . 360 

0.177 

0.016 

32.9 

5.0 

4359. 

1039.44 

0.51 

1051 .07 

12.070 

0.184 

0.01  5 

32.8 

5.0 

4363. 

1040.58 

0.51 

1052.22 

12.084 

0.184 

0.015 

32.6 

5.0 

4580. 

1 102.90 

0.51 

1115.26 

12.810 

0.189 

0.015 

32.4 

5.0 

■4735. 

1  1  48 . 1  5 

0.52 

1161.04 

1 3. 332 

0.192 

0.014 

32.2 

5.0 

4755. 

1  1  54.04 

0.52 

1 1 66.99 

1 3.400 

0.192 

0.014 

32.1 

5.0 

4895. 

1  195.53 

0.52 

1 208 .96 

13.376 

0.195 

0.014 

32.0 

4.0 

3944. 

1 205.98 

0.52 

1219.54 

1  3.995 

0.219 

0.016 

31  .  6 

5.0 

5055. 

1 243. 60 

0.53 

1257.59 

14.413 

0.195 

0.014 

31  .5 

5.0 

5353. 

1 335.01 

0.53 

1 350.07 

15.439 

0.206 

0.013 

31  .  3 

5.  1 

531  6. 

1291 . 38 

0.53 

1 305.92 

14.947 

0.199 

0.013 

31.2 

5.0 

5271  . 

1 309.61 

0.53 

1 324. 37 

15.153 

0.203 

0.013 

31  .0 

5.0 

5  406. 

1 351 .54 

0.53 

1 366. 78 

1 5.627 

0.208 

0.013 

30.8 

5.0 

5  701  . 

1 445.00 

0.54 

1461.34 

1 6. 701 

0.221 

0.013 

30.7 

5.0 

5  707. 

1 446.93 

0.54 

1 463.29 

16.723 

0.221 

0.013 

30.5 

5.0 

5901  . 

1 509.86 

0.54 

1 526.95 

17.451 

0.227 

0.013 

30.4 

4.0 

4735. 

1515.67 

0.54 

1  5  32.8  3 

17.518 

0.255 

0.015 

30.0 

5.0 

6  338. 

1 655.92 

0.55 

1674.  71 

19.096 

0.243 

0.01  3 

29.9 

5.0 

6418. 

1 683. 33 

0.55 

1 702.44 

19.421 

0.253 

0.013 

29.7 

5.0 

6332. 

1 653.88 

0.55 

1 672. 64 

1 9.072 

0.243 

0.013 

29.6 

5.0 

6350. 

1 660.02 

0.55 

1 678.85 

19.144 

0.245 

0.013 

29.5 

5.0 

6462. 

1 698.50 

0.55 

1717.78 

19.605 

0.258 

0.01  3 

29.3 

5.0 

6764. 

1804. 38 

0.55 

1824.89 

20.929 

0.296 

0.014 

29.1 

5.0 

6708. 

1784.51 

0.56 

1  804.79 

20 • 666 

0.28  3 

0.014 

29.0 

5.0 

6757. 

1801  .89 

0.56 

1 822. 37 

20.896 

0.294 

0.014 

28.9 

4.0 

5757. 

1961  .53 

0.5  6 

1 983.87 

23. 378 

0. 444 

0.019 

28.6 

5.0 

6895. 

1851.29 

0.56 

1  8  72.  34 

21 . 589 

0.327 

0.015 

28.5 

5. 1 

7  38  3. 

1977.16 

0.56 

1 999. 68 

23.648 

0 . 400 

0.017 

28.3 

5.0 

7335. 

201 3. 44 

0.57 

20  36.  38 

24. 303 

0.439 

0.018 

28.2 

5.0 

7454. 

2058 . 54 

0.57 

2082.01 

25. 1 38 

0.437 

0.01  7 

28.1 

5.0 

7724. 

21 62.94 

0.57 

2187.63 

26.962 

0.433 

0.016 

27.9 

5.0 

7764. 

21  78 . 66 

0.57 

2203.53 

27.242 

0.445 

0.016 

27.7 

5.0 

7901  . 

2232.98 

0. 57 

2258 . 48 

28.263 

0.498 

0.018 

2  7.6 

5.0 

7778. 

2184.17 

0.57 

2209. 1 0 

27. 341 

0.449 

0.01  6 

2  7.4 

4.0 

6502. 

2324.53 

0.58 

2351 .09 

30.125 

0.599 

0.020 

27.2 

5.0 

8424. 

2447.73 

0.58 

2475. 73 

32. 79  4 

0.599 

0.018 

2  7.1 

4.0 

6541  . 

2344.52 

0.58 

2371.  31 

30.542 

0.609 

0.020 

2  7.0 

5.0 

8  401  . 

2438 .03 

0.58 

2465.91 

32.576 

0.598 

0.018 

26.9 

3.0 

8239. 

2370. 40 

0.58 

2397. 49 

31 .09! 

0.558 

0.018 

2  6.7 

5.0 

8719. 

2574.27 

0.59 

2603. 74 

35.629 

0.632 

0.018 

26.6 

5.0 

8675. 

2555.1 3 

0.59 

2584.  38 

35.194 

0.620 

0.018 

26.5 

5.  1 

8509. 

241 3. 33 

0.59 

2440 . 92 

32.020 

0.588 

0.018 

26. 3 

5.0 

8914. 

2660.  1  7 

0.59 

2690.64 

37.662 

0.671 

0.018 

26.2 

5.0 

8892. 

2650.39 

0.59 

2680. 74 

37.429 

0.669 

0.018 

1  7 


The  array  of  118  altitudes  z.,  measured  densities  p.  and  frac- 

1  i 

tional  statistical  uncertainties  A.  given  in  Table  2.  4  are  the  input  data 
for  the  least  square  mathematical  procedure  described  in  Section  2.  2.  2. 

2.  2.  4  Results  of  Analytical  Fit  for  April  23,  1976  Low  Background 
Betasonde.  Flight 

By  use  of  the  data  in  Table  2.  4  for  z.,  p.  and  A.,  the  least  square 

procedure  was  applied  for  j  =  2  -  5  to  obtain  \he  values  for  the  coeffi- 

m 

cients  A..  For  each  set  of  solutions,  the  results  for  rms  sum  of  the 
weighted?  deviations  between  the  analytical  fit  and  experimental  points  A 
was  determined.  This  is  the  quantity  that  was  minimized  to  produce  the 
solution  obtained.  Ac,  the  rms  unweighted  deviation  sum,  was  also  cal¬ 
culated  from  (2.  14);  both  A  and  Ac  are  to  be  compared  to  Am,  the  rms 
fractional  uncertainty,  from  (2.  12).  Results  for  all  of  these  parameters, 

along,  with  p  from  (2.  17),  are  given  in  Table  2.  5. 
o 

The  rms  fractional  uncertainty,  Am,  is  independent  of  i  and  has 

.  *  m  . 

the  value  0.  039.  As  noted  above,  a  satisfactory  fit  is  one  in  which  the 
difference  between  Ac  and  Am  is  acceptably  small. 


Model  fitting  often  requires  use  of  1  0  t  o  1  4  coefficients  in  a  similar 
procedure  to  obtain  what  is  considered  to  be  an  adequate  fit.  Here,  we 
found  that  for  the  particular  experimental  data  used,  4  or  5  are  suffi¬ 
cient  to  provide  a  fit  whose  accuracy  is  consistent  with  the  accuracy  of 
the  data.  Thus,  as  seen  from  the  table,  A  decreases  rapidly  as  j 

goes  from  2  to  3,  but  with  i  as  large  as  5,  the  accuracy  of  the  lit  is 

111 

marginally  better.  Ac  appears  to  go  through  a  minimum  at  j  =  4.  As 
observed  above,  a  fit  derived  by  minimizing  A,  the  rms  sum  of  the  de¬ 
viations  weighted  inversely  with  the  uncertainties,  is  taken  here  as  the 
proper  procedure  to  follow.  This  causes  those  points  with  most  accuracy 
to  produce  the  most  effect  on  the  fit.  However,  such  a  fit  should  not  be 
allowed  to  produce  a  large  value  of  Ac,  the  rms  sum  of  the  unweighted 
deviations.  In  view  of  the  fact  that  there  is  undoubtedly  some  short-term 
atmospheric  variability  in  density,  it  is  unlikely  that  a  fit  in  which  Ac 
approaches  Am  much  more  closely  than  in  the  present  case  would  be 
expected.  Thus,  in  the  present  instance  at  least,  a  criterion  for  select¬ 
ing  the  most  satisfactory  of  the  fits  obtained  by  minimizing  A,  is  that 

particular  one  for  which  Ac  is  a  minimum;  i  =  4  here.  Whether  such 

■  m 

a  dependence  on  j  will  always  occur  is  not  presently  known.  When  the 
results  of  additional  flights  are  available,  futher  investigation  of  t  he 
method  of  selecting  the  optimum  fit  should  he  made. 


Table  2.  5  Parametric  Results  of  Analytical  Fit 
for  LBB  Flight  on  April  23,  1976 


The  quantity  FT(z )  given  by  (2.  18)  is  actually  !v>,he  used  in  (2.  ]  0)  to  pro¬ 
vide  the  optimum  analytical  fit  to  the  data.  However,  for  purposes  of  com¬ 
parison  with  models  it  is  useful  to  have  the  value  of  H(zl.  This  tan  he  de¬ 
termined  as  follows.  Note  from  (2.  9)  that 


Use  of  (2.  IS)  then  yields 

lit/) 


■  m 


y  E  a.o  - 1  )(v.-v.ny 


_  7 


(2.  1  9) 


(2.20) 


(2.21) 


Consequently,  once  the  Aj  have  been  determined  it  is  possible  to  calc  ulate 
tlie  value  of  H(z)  as  well  as  ll(z). 

For  the  fit  with  j  -  4,  results  are  given  in  Table  2.  6  and  Figures 
2.2  and  2.  3.  The  first  three  columns  of  the  table  give  z-,  f,-  and  (ALT, 
DENSITY  and  UNC  in  Table  2.4).  The  next  three  give  p(z),  H(z),  and  H(z) 
as  calculated  from  (2.10),  (2.18),  and  (2.21),  respectively.  The  nc-xt  column 
gives  the  fractional  difference  between  measured  pj  and  calculated  p(z). 

The  final  column,  labeled  RMS  FRAC,  is  a  running  rms  average  of  the 
fractional  difference  for  that  altitude  with  that  of  the  four  previous  (higher 
altitude  points).  Hence,  it  is  assigned  the  value  zero  for  the  four  highest 
altitude  points. 


Examination  of  the  variation  of  RMS  FRAC  with  altitude  shows  those 
regions  in  which  the'  best  fit  is  obtained.  At  the'  highest  altitudes  the  aver¬ 
age  difference  for  any  single  point  is  as  high  as  7-S"i.,  due  to  the  small 
number  of  counts  recorded  there.  As  the  sonde  descends  il  goes  through 
a  region  where  the  average  difference-  is  on  the  order  of  2"u.  There  is, 
however,  a  lower  altitude  region  (below  about  20  km)  in  which  beta  par¬ 
ticle  energy  loss  dominates  the  effects  of  scattering,  and  the  count  rat  o 
deer  ea  s  <-s  .is  z  decreases  (Ref,  1.  !)•  Thus,  iiea  r  20  km  there  is  little 
variation  of  count  rate  with  z,  and  the  density  uncertainty  becomes  very 
large  in  spite  of  the  large  count  rates.  This  region  is  being  approached 
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Table  2.6  Altitude  Dependence  Results  of 
Analytical  Fit  for  LRU  Flight  on 
April  2  L  1  '76 


j 


m 


! 


results  of  fit 

J  A(J) 

1  -0. 35753E+01 

2  0.1701 7E+00 

3  -0.76747E-03 

4  -0 . 1 8  776E-05 

Z0-  26.2  KM 
FHO0=  35.707  G/M**  3 


RMS  FRAC=RMS  AVG  OF  PREV  5  DIFF  FRACS 


ALT 

RHO-MEAS 

UNCER 

RHO-CALC 

HBARCZ ) 

H  ( Z  > 

DIFF 

KM 

G/M**  3 

G/M**  3 

G/M**  3 

KM 

KM 

FRAC 

64.0 

0.168 

0.019 

0.190 

7.221 

9.606 

-0.117 

64.3 

0.198 

0.021 

0.184 

7.236 

9.661 

0.073 

64.2 

0.182 

0.020 

0.186 

7.231 

9.642 

-0.024 

63.5 

0.185 

0.020 

0.201 

7.198 

9.516 

-0.077 

62.9 

0.232 

0.023 

0.214 

7.170 

9.412 

0.086 

62.1 

0.225 

0.023 

0.233 

7.1  33 

9.276 

-0.033 

61  .0 

0.291 

0.026 

0.262 

7.083 

9.096 

0.109 

60  *  0 

0.298 

0.026 

0.29  3 

7.0  38 

8.940 

0.017 

59.1 

0.357 

0.029 

0.324 

6.998 

8.805 

0.100 

57.9 

0.355 

0.028 

0.372 

6.947 

8.632 

-  0 . 046 

55.9 

0.504 

0.034 

0.471 

6.862 

8.360 

0  •  070 

54.9 

0.553 

0.036 

0.532 

6.821 

8.232 

0 . 040 

53.9 

0.597 

0.037 

0.601 

6.781 

8.108 

-0 . 006 

53.0 

0 . 644 

0.0  38 

0.672 

6.745 

8.001 

-0.041 

52.0 

0.839 

0.044 

0.762 

6.706 

7.885 

0.101 

51  .2 

0.854 

0.044 

0.844 

6.675 

7.796 

0.012 

50.5 

0.889 

0 .046 

0.923 

6.648 

7.719 

-0.037 

50.0 

0.897 

0.051 

0.985 

6.629 

7.666 

-0 . 090 

49.5 

1  .069 

0.056 

1 .052 

6.611 

7.61  3 

0.016 

48.6 

1  .290 

0.055 

1.185 

6.577 

7.521 

0.089 

48.1 

1.191 

0.053 

1 .267 

6.559 

7.471 

-0.060 

47.7 

1 .316 

0.055 

1 . 337 

6.544 

7.431 

-0.015 

47.2. 

1 .388 

0.064 

1 .430 

6.526 

7.38  3 

-0.029 

46.8 

1.443 

0.065 

1.510 

6.512 

7.344 

-0.044 

46.4 

1  .437 

0 . 064 

1 .594 

6.498 

7.306 

-0.099 

45.3 

1  .835 

0.104 

1 .731 

6.476 

7.250 

0 . 060 

45.3 

1 .945 

0.075 

1  .855 

6.459 

7.205 

0.048 

45.0 

1  .871 

0.074 

1 .935 

6.448 

7.178 

-0.033 

44.2 

2.044 

0.075 

2.164 

6.421 

7.  106 

-0.055 

43.9 

2.239 

0.0  79 

2.257 

6.410 

7.080 

-0 . 008 

43.4 

2.510 

0.076 

2.423 

6.393 

7.037 

0.036 

43.  1 

2.322 

0.0  72 

2.529 

6.  38  3 

7.012 

-0.082 
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RMS 
FRAC 
0.000 
0.000 
0 . 000 
0.000 
0.081 
0.0  64 
0.073 
0.073 
0.078 
0.071 
0.077 
0.062 
0.061 
0.046 
0.061 
0.052 
0.052 
0.066 
0.063 
0.060 
0.065 
0.063 
0.051 
0.054 
0.057 
0.057 
0.061 
0.061 
0.063 
0.045 
0.040 
0.049 


Table  .  (>  (cont'd) 


ALT 

RHO-MEAS 

KM 

G/M** 3 

42.8 

2.622 

42.5 

2.699 

42.2 

2.821 

41.9 

2.831 

41.7 

3.086 

41 . 1 

3.410 

40.8 

3.499 

40.5 

3.671 

40.2 

3.666 

40.0 

3.931 

39.7 

4.042 

39.5 

4.285 

39.3 

4.451 

39.0 

4.671 

38.5 

4.768 

38.3 

5.008 

38. 1 

5.225 

37.9 

5.618 

37.6 

5.654 

37.4 

5.916 

37.2 

6.054 

37.0 

6.096 

36.8 

6.498 

36.2 

7.242 

36.  1 

7.693 

35.9 

7.725 

35.7 

7.478 

35.5 

8.097 

35.4 

8.086 

35.2 

8.286 

35-  1 

8.819 

34.9 

8.743 

34.5 

9.226 

34-4 

9.692 

34.3 

9.388 

34. 1 

9.123 

34.0 

9.811 

33.9 

10.267 

33.7. 

10.332 

33.6 

11.012 

33.5 

1 1 .227 

33.1 

1  1.205 

33.0 

11.360 

32.9 

12.070 

52.8 

12.084 

32.6 

12.810 

UNCER  RHO-CALC 
G/M**  3  G/M**  3 

0.078  2.639 

0.079  2.755 

0.080  2.877 

0.091  3.004 

0.095  3.093 

0.091  3.375 

0.092  3.526 

0.095  3.685 

0.094  3.852 

0.098  3.967 

0.103  4.148 

0.106  4.273 

0.108  4.402 

0.123  4.604 

0.123  4.962 

0.114  5.114 

0.116  5.271 

0.120  5.433 

0.120  5.685 

0.123  5.861 

0.122  6.043 

0.123  6.230 

0.128  6.424 

0.136  7.045 

0.141  7.154 

0.141  7.379 

0.139  7.611 

0.145  7.851 

0.145  7.974 

0.147  8.226 

0.153  8.356 

0.170  8.621 

0.156  9.179 

0.163  9.324 

0.159  9.472 

0.156  9.775 

0.164  9.930 

0.171  10.088 

0. 172  10.41 3 

0.177  10.579 

0.199  10.748 

0.178  11.454 

0.177  11 .638 

0.184  11.825 

0.184  12.015 

0.189  12.406 


HBAR(Z  >  H ( Z ) 
KM  KM 

6.373  6.986 

6.363  6.961 

6.353  6.936 

6.343  6.912 

6.336  6.895 

6.316  6.847 

6.307  6.823 

6.297  6.800 

6.287  6.776 

6.281  6.761 

6.271  6.738 

6.264  6.722 

6.258  6.707 

6.248  6.684 

6.233  6*647 

6.226  6.632 

6.220  6.618 
6.214  6.603 

6.204  6.582 

6.198  6.567 

6.192  6.553 

6.186  6.539 

6.180  6.525 

6.161  6.483 

6.158  6.476 

6.152  6.462 

6.146  6.448 

6.140  6*435 

6.137  6.428 

6.131  6.414 

6.128  6.408 

6.122  6.394 

6.110  6.368 

6.107  6.361 

6.104  6.355 

6.098  6.341 

6.095  6.335 

6.092  6.328 

6.086  6.315 

6.083  6.309 

6.080  6.303 

6.068  6.277 

6.066  6.271 

6.063  6.264 

6.060  6.258 

6.054  6.245 


DIFF  PM  $ 

FRAC  FRAC 

-0.007  0.047 

-0.02!  0.041 

-0.019  0.042 

-0.058  0.047 

-0.002  0.029 

0.010  0.029 

-0.008  0.028 
-0.004  0.027 

-0.048  0.022 

-0.009  0.023 

-0.025  0.025 

0.003  0.025 

0.011  0.025 

0.015  0.015 

-0.039  0.02? 

-0.021  0.021 

-e.009  0.022 

0.034  0.026 

-0.006  0.025 

0 . 009  0.019 

0.002  0.017 

-0.022  0.019 

0.012  0.012 

0.028  0.017 

0.075  0.038 

0 . 047  0.043 

-0.017  0.043 

0.031  0.045 

0.014  0.043 

0.007  0.027 

0.055  0.030 

0.014  0.020 

0.005  0.027 

0.039  0.031 

-0.009  0.031 

-0.067  0.036 

-0.012  0.035 

0.018  0.036 

-0.008  0.032 

0.041  0.036 

0.045  0.029 

-0.022  0.030 

-0.024  0.031 

0.021  0.032 

0.006  0.026 
0.033  0.023 


Table  2.  6  (cont'd) 


ALT 

RHO-KEAS 

UNCER 

RHO-CALC 

HEAP (  Z  ) 

HC2) 

DIFF 

RMS 

KM 

G/M**3 

G/M**  3 

G/M** 3 

KM 

KM 

FRAC 

FRAC 

32.4 

1 3.332 

0.192 

12.810 

6.048 

6.233 

0.041 

0.027 

32.2 

1 3.400 

0.192 

1 3.228 

6.042 

6.221 

0.013 

0.026 

32.  1 

1 3.876 

0.195 

1 3.443 

6.039 

6.214 

0.032 

0.028 

32.0 

1 3.995 

0.219 

1 3. 661 

6.037 

6.208 

0.024 

0.0  30 

31.6 

14.413 

0.195 

14.572 

6.025 

6.184 

-0.011 

0.027 

31.5 

15.439 

0.206 

14.810 

6.022 

6.178 

0.043 

0.027 

31 . 3 

14.947 

0.199 

15.297 

6.017 

6.1  65 

-0.023 

0.029 

31.2 

15*153 

0.20  3 

15.548 

6.014 

6.159 

-0.025 

0.027 

31 .0 

15.627 

0.208 

16.061 

6.008 

6.147 

-0.027 

0.028 

30.8 

16.701 

0.221 

16.593 

6.002 

6.  1  35 

0.007 

0.027 

30.7 

16.723 

0.221 

16.866 

6.000 

6.129 

-0.008 

0.020 

30.5 

1 7.451 

0.227 

1 7.426 

5.994 

6.117 

0.001 

0.017 

30.4 

17.518 

0.255 

17.713 

5.991 

6.112 

-0.011 

0.014 

30.0 

19.096 

0.243 

18.914 

5.980 

6.088 

0.010 

0.008 

29.9 

19.421 

0.253 

19.227 

5.977 

6.082 

0.010 

0.009 

29.7 

19.072 

0.243 

19.870 

5.972 

6.071 

-0.040 

0.020 

29.6 

19.144 

0.245 

20.200 

5.969 

6.065 

-0.052 

0.031 

29.5 

19.605 

0.258 

20.536 

5.966 

6.059 

-0.045 

0.036 

29.3 

20.929 

0.296 

21 .226 

5.960 

6.047 

-0.014 

0.037 

29.  1 

20.666 

0.283 

21.941 

5.955 

6.036 

-0.058 

0.045 

29.0 

20.896 

0.294 

22.307 

5.952 

6.030 

-0.063 

0.050 

28.9 

23.378 

0.444 

22.681 

5.949 

6.025 

0.031 

0.046 

28.6 

21.589 

0.327 

23.840 

5.941 

6.008 

-0.094 

0.059 

28.5 

23-648 

0.400 

24.241 

5.9  38 

6.002 

-0.024 

0.060 

28.3 

24.303 

0.4  39 

25.063 

5.933 

5.991 

-0.0  30 

0.055 

28.2 

25.138 

0.437 

25.485 

5.9  30 

5.985 

-0.014 

0.048 

28. 1 

26.962 

0.433 

25.914 

5.927 

5.980 

0 . 040 

0.050 

27.9 

27.242 

0 . 445 

26.797 

5.922 

5.969 

0.017 

0.027 

27.7 

28.263 

0.498 

27.71 1 

5.917 

5.957 

0.020 

0.026 

27.6 

27.341 

0 . 449 

28.180 

5.914 

5.952 

-0.0  30 

0.026 

27.4 

30.125 

0.599 

29.144 

5.909 

5.941 

0.0  34 

0.029 

27.2 

32.794 

0.599 

30.143 

5.903 

5.930 

0.088 

0.046 

27.  1 

30.542 

0.609 

30.655 

5.900 

5.925 

-0.004 

0.045 

27.0 

32.576 

0.598 

31.178 

5.898 

5.919 

0.045 

0.049 

26.9 

31 .091 

0.558 

31 . 709 

5.895 

5.914 

-0.019 

0 . 047 

26.7 

35.629 

0.632 

32.801 

5.890 

5.903 

0.086 

0.059 

26.6 

35.194 

0.620 

33.361 

5.887 

5.898 

0.055 

0.051 

26.5 

32.020 

0.588 

33.932 

5.884 

5.892 

-0.056 

0.057 

26.3 

37.662 

0.671 

35.105 

5.879 

5.882 

0.073 

0.062 

26.2 

37.429 

0.669 

35.707 

5.876 

5.876 

0 .048 

0.065 

TOT.  RMS  FRAC.  RHOMEAS-RHOCALC  DI FF-  0.04475 
TOT.  RMS  FRAC.  RHOMEAS.  UNCERTAINTY*  0.03909 


Altitude  (km) 

50  5  5 


(.0 


Altitude  (kill) 

.  5  l.nw  Bae  ky;  r  mind  Belnsmide  Air  Ihnsity  Result:-  1 1  >  r  l-'li.'.lil 
.  1 1  White  Sands  Missile  K.iime  mi  -1  >  /  V  t  ■ .  -1  •  -  « >  •  i  !  m. 


at  the  lowest  altitudes  in  Table  2.6,  which  accounts  for  much  of  the  aver¬ 
age  difference  of  about  5%  there.  As  noted  previously,  the  4  measure¬ 
ments  made  in  the  vicinity  of  the  64  km  apogee  during  a  period  of  about 
20  seconds  are  particularly  important,  since  the  sonde  descends  rapidly 
for  some  distance  ther  eafter.  The  data  from  these  four  points  could  have 
been  averaged  together  to  give  one  measurement  at  (he  average  altitude  of 
64.0  km.  The  result  would  have  been  <  > .  1  S  -i  -i  g/m^  with  a  fractional  uncer¬ 
tainty  of  .054  (about  half  that  of  each  of  the  four  points).  i'he  calculated 
result  from  the  fit  is  .190  g/m  ,  a  fractional  difference  of  only  .036. 

Thus,  these  four  points,  as  utilised  with  least  square  procedure,  serve 
to  define  p(x)  accurately  in  the  apogee-  altitude  region. 

It  is  of  interest  to  compare  the  ll(z)  values  with  those  of  an  approp¬ 
riate  model  atmosphere,  such  as  that  of  (lie-  Spr  ing/Fall  Mid- Latitude-  Model 
(Ref.  2.4).  Values  of  II(z)  are  not  given  in  that  particular  model.  It  is 
straightforward,  however,  to  derive  the  associated  !!(/.)  by  fitting  the  two 
p(y.)  values  that  bracket  the  altitude  of  interest.  13y  use  of  this  procedure 
and  the  results  of  Table  2.6  we  obtain  tlu-  following: 


Table  2.  7 

Comparison  of  Scale  Heights 
for  4/23/76  l-'light 
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-  - 
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4  0 
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4  5 
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i  .  IS 

50 

K.  0-1 

7.  i-7 

S 

5.  40 

s.  .  t 

t)0 

h.  1 1 

S.  0.1 

ti^ 

H .  OS 

~o,  70 

In  general  the  ll(z)  values  from  the  analytical  lit  arc  .-  lightly  below  those 
ot  the  model  in  the  2s- ’-t'  km  region,  close  in  the  3s-1-'-  1  in  region,  and 
somewhat  higher  in  the  Ml-(>6  km  region.  Certainly,  the  agreement  is 
i  r\  reasonable  on  the  average. 
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2.  2.  5  Discussion 

In  summary,  we  believe  that  the  results  oi  applying  the  least  square 
mathematical  procedure  developed  in  Section  2.  2.  2  to  the  LBB  flight  data  in 
Section  2.2.3  show  the  procedure  to  be  valid.  An  accurate  analytical  lit  was 
obtained  that  is  capable  of  representing  the  entire  tabular  array  of  data, 
and  we  believe  that  the  procedure  should  be.  applied  to  future  Lilli  or  lieta- 
sonde  11  flights. 

It  should  be  observed  that  a  similar  procedure  can  also  be  applied  to 
other  altitude-dependent  measurements,  for  example  temperature  and  pres¬ 
sure,  provided  a  tabular  array  of  the  measured  quantity,  and  its  uncertainty, 
is  available.  In  the  case  of  pressure  measurements  the  above  procedure 
could  be  used  almost  in  its  exact  form,  since  pressure  also  depends  expo¬ 
nentially'  on  altitude.  For  temperature  the  analytical  function  would  have 
to  be  chosen  to  represent,  in  an  average  manner,  the  known  variation  of 
that  quantity  with  altitude. 
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3.  BETASONDE  II  DESIGN  AND  OPERATION 


In  this  section  the  design  of  the  advanced  model,  Betasonde  II,  is  given, 
and  details  of  the  unit,  as  fabricated,  are  presented.  Betasonde  II  is  present¬ 
ly  available  for  flight  and  could  readily  be  fitted  into  rockets  other  than  the 
Areas. 

3.  1  Description 

A  photograph  of  the  completed  sonde,  which  was  originally  designed  for 
flight  on  an  Areas  rocket,  is  shown  in  Fig.  3.1.  The  basic  configuration  is  given 
in  Fig.  3.2.  A  tripod-supported  boom  carries  the  annular  beta  source  which  ir¬ 
radiates  the  air  volume  in  the  region  between  source  and  detector,  outside  the 
shield.  The  forward-scattered  beta  rays  are  detected  by  a  surface  barrier 
type  semiconductor  detector  located  at  the  base  of  the  tripod.  A  shield  is 
placed  between  the  source  and  detector  to  prevent  direct  passage  of  beta  par¬ 
ticles  from  the  source  to  the  detector.  A  plastic  disc  is  placed  on  the  side  of 
the  shield  that  faces  the  source,  in  order  to  reduce  hr ems strahlung  generation. 
The  distances  of  the  source  to  the  detector  and  to  the  shield  are  fixed  for  opti¬ 
mum  system  response  characteristics.  Work  to  optimize  that  response  was 
carried  out  during  the  development:  phase  of  this  contract  and  is  contained  in 
our  report:  Atmospheric  Density  Measurement  In  The  Middle  Atmosphere 
(Ref.  1.1  ). 

The  original  concept  of  the  high-altitude  Betasonde  utilized  an  annular 
surface  barrier  detector  with  the  source  boom  mounted  through  the  center 
of  that  detector.  Although  this  arrangement  has  the  advantage  of  simpler 
mechanical  construction,  it  has  been  found  that  the  cost  of  a  rnggedized, 
light  weight,  annular  detector  outweighs  these  advantages.  The  choice  of 
the  below  described  non-annular  detector  is  based  on  its  cost  effectiveness 
and  on  the  manufacturer's  experience  with  the  construction  of  such  a  special 
dot  eel  or. 


Fig.  3.1  Photograph  of  Completed 
Betasonde  II  Atmospheric 
Density  Sonde 
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The  detector  is  a  ruggedized  surface  barrier  type  with  an  active  area 
of  450  mm^,  sufficient  to  give  a  reasonable  count  rate  up  to  an  altitude  of 
80  km.  Since  a  standard  detector  inherently  responds  to  light  in  the  visible 
spectrum,  it  is  necessary  for  this  application  to  have  a  specially  processed 
detector  featuring  a  light-tight  front  surface.  According  to  the  manufacturer 
(Ortec,  Inc.)  a  light-tight  front  surface  can  be  obtained  with  120pg/cm^  of 
aluminum,  which  is  vapor  deposted  directly  on  the  detector  material.  The 
latter  is  p-type  silicon  (rather  than  the  standard  n-type)  in  order  for  the 
aluminum  to  become  the  rectifying  contact  of  the  surface  barrier  diode,  and 
thus  the  front  of  the  detector.  Other  features  of  the  detector  arc  low  noise 
(18  keV  max  for  beta  particles)  and  a  depletion  depth  of  300pm  at  a  bias  of 
100V.  The  relatively  low  bias  voltage  is  a  desirable  feature,  simplifying 
the  design  and  assuring  low  power  dissipation  of  the  bias  voltage  supply. 

The  source  consists  of  an  annular  disc  with  the  radioisotope  (Pm- 147) 
deposited  and  sealed  on  the  underside  of  it.  The  center  hole  of  the  disc  is 
tapped  for  convenient  mounting  to  the  boom.  A  cone  fabricated  from  alum¬ 
inum  serves  as  a  protective  cover  and  firmly  holds  the  source  in  place. 

Source  intensity  is  of  the  order  of  ^00  mCi  for  .SO  km  maximum  altitude. 

The  lower  portion  of  the  sonde  contains  all  necessary  electronics  to 
process  the  signals  from  the  detector  and  to  supply  the  required  power  for 
the  electronics  as  well  as  the  detector.  The  system  circuitry  is  housed  in 
a  cylindrical  container,  constructed  of  aluminum  with  a  diameter  of  3-1/2  in. 
by  6-3/8  in.  long.  The  charge  sensitive  preamplifier  (CSPA)  is  additionally 
shielded  by  its  own  housing  to  protect  it  from  electro-magnetic  interference. 
The  rest  of  the  electronics  is  mounted  on  three  pr inter!  circuit  hoards  which 
individually  plug  into  a  motherboard.  The  internal  power  supplies  (DC/ DC 
converters)  are  mounted  to  the  bottom  cover  for  improved  heat  transfer  to 
the  outside  environment. 

The  telemetry  antenna,  not  shown  in  Fig.  3.1,  may  be-  attached  to  the  boom 
and  along  one  leg  of  the  tripod.  Grommets  in  the  top  and  bottom  of  the  elec¬ 
tronics  housing  are  provided  for  a  coaxial  antenna  feed  line  which  would  pass 
through  the  housing.  The  mechanical  interface  of  the  rocket  telemetry  sec¬ 
tion  and  the  sonde  simply  consists  of  four  threaded  spacers  used  to  holt  the 
two  units  together.  Electrically  the  two  units  are  hard-wired  together. 

3.  2  '  Opcrat ion 

The  electronic  system  of  the  sonde  is  illustrated  in  block  diagram  form 
in  Fig.  3.3.  The  instrument  consists  essentially  of  a  solid  state  det  eet  or /amp¬ 
lifier/threshold  discriminator  chain  driving  selectively  a  binary  divider  or 
passing  right  through  to  be  applied  to  telemetry  (I'M)  after  suitable  shaping 
and  multiplexing  with  detector  temperature  information. 


tate  Detector,  light  tight 


The  two  signals  (pulses  or  temperature)  are  alternately  applied  to  TM 
by  means  of  the  multiplexer.  A  more  detailed  description  of  the  system  is 
given  below. 

The  forward-scattered  beta  rays  incident  on  the  solid  state  detector  pro¬ 
duce  hole-electron  pairs  in  the  silicon  detector  (1  hole-elect  ran  for  each  2.3 
eV  of  deposited  energy)  which  are  swept  up  by  the  electric  field  created  by  the 
bias  voltage,  resulting  in  a  charge  pulse  proportional  to  the  incident  deposited 
energy  at  the  detector  output.  Overbiasing  the  solid  state  detector  to  t  1  10Y 
assures  that  it  is  totally  depleted  and  will  remain  so  under  varying  ambient 
conditions;  and  all  charge  will  be  collected. 

The  charge  pulse  is  converted  to  a  voltage  pulse  by  the  charge  sensitive 
preamplifier  (CSPA)  and  is  then  amplified  and  shaped  by  the  shaping  amplifiers. 
Double  differentiation/ double  integration  shaping  is  used  for  near  optimum  sig¬ 
nal  to  noise  ratio.  Use  of  a  charge  sensitivie  preamplifier  assures  that  the 
shaping  amplifier  output  conversion  gain  of  4  volts  per  MeV  of  deposited  ener¬ 
gy  will  be  insensitive  to  variation  in  solid  slate  detector  or  stray  capacitances. 

An  output  pulse  is  produced  by  the  threshold  discriminator  (Th.  Disc.)  each 
time  the  shaping  output  level  exceeds  the  fixed  reference  level.  This  reference 
is  set  at  a  voltage  equivalent  to  about  2.5  times  the  energy  noise  level  of  the  de¬ 
tector.  As  stated  above,  the  maximum  FWI1M  noise  output  from  the  detector 
is  18  keV  -  including  electronics  noise  -  which  results  in  an  energy  threshold 
level  of  45  keV.  The  discriminator  pulses  are  standardized  as  to  pulse  width 
and  height  with  a  monostable  multivibrator  or  One-Shot  (O/S). 

The  One-Shot  output  pulses  are  simultaneously  applied  to  a  binary  Divider 
(•*■  8),  an  OR  gate  (Select)  and  a  Digital  Rut  cinder.  As  long  as  the  O/S  pulse 
rate  does  not  exceed  about  10  keps  the  O/S  pulses  are  directly  mu  put  to  the 
TM  via  the  OR  gate  to  give  the  desired  air  density  information.  If  (lie  count 
rate  exceeds  about  10  keps,  the  output  is  selected  from  the  -MS  circuit,  thus 
extending  the  maximum  possible  input  count  ing  rate  by  a  factor  of  8.  The 
selection  process  is  based  on  the  ratemeter  output  and  its  associated  control 
and  logic  circuits.  The  state  of  the  ratemeter  out  put  is  continuously  sampled. 
When  the  input  count  rate  exceeds  10,240  ops  the  output  goes  low  and  the  con¬ 
trol  circuit  selects  the  -j-8  output  as  the  air  density  information  channel  for 
TM.  Relow  that  count  rate  threshold  the  O/S  pulses  are  directly  routed  to 
TM  via  the  Select  gate.  The  Clock  and  Delay  circuits  associated  with  the 
Ratemeter  and  Control  provide  the  necessary  logic  signals  for  the  sampling 
operation  and  furnish  a  time  base  for  the  Ratemeter. 

The  temperature  of  the  solid  state  detector  is  continuously  measured 
with  a  small  thermistor  bead  in  contact  with  the  detector  housing.  The  out¬ 
put  voltage  of  the  thermistor  varies  according  to  temperature  in  a  non-linear 
fashion  and  is  processed  by  a  voltage  to  frequency  converter  in  the  Tempera¬ 
ture  Monitor  circuit.  Full  scale  output  is  200  llz  at  4  r-> 0 ‘ ’ C '  non-linearly 
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decreasing  to  about  40  IIz  at  -10°C.  This  range  gives  ample  margin  for  the 
expected  temperature  extremes  of  0°  to  +  45°C  derived  from  measurements 
of  other  payloads  flown  on  Areas  rockets.  The  range  also  covers  the  upper 
operating  temperature  limit  of  the  detector  of  45°C. 

Since  the  TM  transmitter  can  only  handle  one  information  channel,  the 
selected  pulse  output  and  the  Temperature  Monitor  output  are  multiplexed. 

Time  allocation  is:  50  seconds  for  pulse  output,  5  seconds  for  temperature 
output.  The  FM/FM  transmitter  can  process  signal  rates  of  up  to  10  kHz. 

To  achieve  sufficient  bandwidth  for  a  pulse  output,  a  dead  time  circuit  in  the 
Pulse  Shaper  limits  the  pulse  repetition  rate  to  5  kHz.  Since  the  sonde  input 
pulse  rate  can  exceed  that  rate,  up  to  10  kHz,  a  dead  time  correction  must 
be  applied  as  described  in  the  annual  report  (Ref.  1.1).  Also  included  in  the  Pulse 
Shaper  is  a  level  translator  converting  standard  logic  pulses  to  a  0  to  -12V 
pulse  of  40/US  duration.  This  makes  the  output  compatible  with  FM/FM  trans¬ 
mitters  previously  flown  on  Areas  rockets. 

All  necessary  power  for  the  system  is  derived  from  2  DC/ DC  converters 
connected  to  the  +28V  battery  pack  of  the  TM  section.  The  low  voltage  con¬ 
verter  delivers  _+12V  at  +22  ma  and  -16  ma,  respectively.  The  H.  V.  conver¬ 
ter  delivers  +11  0V  at  max  4  ma.  The  power  requirement  for  the  battery  is 
28V  at  120  ma  for  the  above  output  power. 

For  completeness,  the  detailed  schematics  of  t  lie  various  circuits  of  the 
block  diagram  are  included  in  this  report,  and  are  contained  in  Figs.  3.4  thru 
3.7.  These  show  the  as-constructed  circuits  contained  in  the  completed  sonde 
in  Figs.  3.1  thru  3.3. 

3.  3  Prometheum  Radiation  Source 

A  Pm-147  source  of  intensity  500  mCi  was  ordered  from  Isotope  Prod¬ 
ucts  Labs  in  October,  1977.  When  received,  however,  our  measurement  s 
indicated  an  actual  emission  less  than  100  mCi.  It  was  returned  to  IPL  for 
possible  repair,  but  when  returned,  it  was  of  the  same  intensity.  Thus, 
all  measurements  reported  above  were  made  with  the  older  Pm-147  source, 
which  had  an  emission  near  100  mCi  when  used.  We  believe  it  likely  that 
the  specific  activity  of  the  isotope  used  to  prepare  the  IPL  source  was  not 
as  high  as  necessary.  We  feel  certain  that  a  source  with  >  500  mCi  Pm-147 
emission  can  be  fabricated.  But  it  will  require  more  procurement  effort 
and  quality  assurance  than  was  possible  during  this  limited  program. 
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B.;i  a  sonde  II  -  Schematic,  Shaping  Amplifier  Thresh 
Discr  iminat  or  ,  One-Shot 


Beiasondc  II  -  Schematic,  Clock,  Delay, 


3.  4  Specifications 


The  pertinent  parameters  of  the  Betasonde  II  atmospheric  density  sond<. 

are  tabulated  below: 

Altitude  Range*  25-65  km 

Beta  Source*  Pm-147,  approx.  100  mCi 

Detector  Surface  barrier  p-type  silicon  detect  or,  light-tight, 

450  mm  ,  300  gm  depletion  depth. 

Telemetry  Output  Two  multiplexed  outputs:  Air  density  pulses  out  fur 

50s,  detector  temp,  pulses  out  for  5s;  4gs  pulses, 

0  to  -12V  of  5  keps  max.  repetition  rate  v>,< 

Size  3-1/2  in.  diameter  by  6-3/8  in.  long  electronics 

housing  with  integral  boom  extending  8-1/2  in. 
above  housing. 

Weight  Total  1.87  lbs  (850  g) 

Power  28V  at  120  ma  from  battery  external  to  betasonde 

*The  altitude  range  can  be  significantly  extended  upward  by  use  of  a  more 
intense  beta  source.  For  a  500  mCi  source  the  altitude  maximum  would 
be  80-90  km  depending  upon  the  trajectory  characteristics. 

**This  implies  that  if  betasonde  is  used  with  other  than  standard  Areas  1680  MU 
transmitter,  the  TM  channel  used  must  be  capable  of  5  kHz  maximum  trans¬ 
mission.  Capability  to  count  the  pulses  is  all  that  is  required,  so  that  some 

pulse  shape  distortion  by  the  transmitter  can  be  tolerat  ed. 


4.  SUMMARY  AND  CONCLUSIONS 


During  Phase  I,  Research  and  Development,  the  following  was  ac¬ 
complished. 

1.  A  technique  was  developed  for  making  an  analytic  fit  to  an 
entire  array  of  density  versus  altitude  data  on  a  least-square- 
basis.  The  procedure  was  applied  to  data  from  a  Low  Back¬ 
ground  Betasonde  flight  with  ,Tood  results. 

2.  By  use  of  theoretical  calculations  and  experimental  measure¬ 
ments,  it  was  shown  that  Pm-147  is  the  optimum  choice  source 
for  the  single- scattering  betasonde.  This  conclusion  would  be 
modified  if  it  were  possible  to  reduce  the  energy  threshold  (that 
energy  above  which  all  betas  are  counted),  to  about  10-20  keV. 

In  that  case,  Ni-63  would  probably  be  optimum,  and  the  maxi¬ 
mum  altitude  could  be  above  100  km. 

3.  A  laboratory  model  of  Betasonde  II,  the  semiconductor -del  ector 
version,  was  constructed  and  calibrated  with  Pm-147.  Although 
some  wall-effect  count  rate  was  evident  at  very  low  densities 
(due,  apparently,  to  the  open  geometry  of  the  detector),  it  was 
possible  to  obtain  meaningful  calibration  data. 

4.  The  calibration  data  were  used  to  show  that  with  a  Pm-147 
source  of  about  500  mCi,  it  should  be  possible  to  obtain  about 
+  6%  accuracy  in  the  90  km  region  for  an  Areas  launch  having 
apogee  near  that  altitude.  Below  that  altitude  the  accuracy 
would  be  better. 

During  Phase  II,  Design  and  Fabrication,  the  following  was  accom¬ 
plished. 

1.  The  design  of  the  electronics  circuitry  was  completed.  This  in¬ 
cludes  the  preamplifier,  analog  processing,  digital  logic,  timing  and  output 
pulse  shaping  circuitry.  The  total  power  requirement  was  determined  as 
given  in  the  Specifications,  Section  3. 

Z.  Breadboards  of  all  circuits  were  constructed  in  order  to  verify 
the  design  requirements  and  goals  over  the  expect  ed  temperature  r.initc. 

3.  The  procurement  of  the  annular  Pm-147  source  from  Isotope 
Products  Labs  was  prepared  and  was  confirmed  on  October  17,  1  977. 

When  the  source  was  delivered,  it  did  not  have  the  500  mCi  emission  ex¬ 
pected.  A  proper  source  can  be  obtained,  but  would  require  more  time 
and  effort  than  could  be  applied  under  the  present  program. 
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4.  During  the  final  portion  of  Phase  II  the  Betasonde  was  fabricated, 
tested  and  calibrated.  It  meets  all  of  the  original  design  objectives  and, 
with  proper  radiation  source,  would  be  capable  of  providing  density  alti¬ 
tude  profiles  to  at  least  80  km  -  90  km. 

As  described  in  Section  ],  the  final  phase  of  this  program  would  be 
Flight  and  Data  Analysis.  During  this  phase  of  the  work  (not  presently 
funded  by  the  Army  Research  Office)  the  instrument  must  be  integrated 
into  an  Areas  (or  other  rocket)  payload  for  launch,  and  calibrated  at  both 
Panametrics  and  in  a  large  (40'  -  60')  chamber  in  order  to  verify  the  high 
altitude  (>  70  km)  portion  of  the  calibration  curve.  The  Betasonde  II 
should  then  be  flown  to  the  80-90  Jem  region  at  least  twice,  and  the  analy¬ 
tical  fitting  procedure  described  above  should  be  applied  to  the  data.  This 
will  form  the  basis  for  routine  direct  measurement  of  atmospheric  density 
up  to  at  least  the  80-90  krn  region,  with  possible  extension  to  100  km  by 
use  of  more  intense  sources  or  lower  electronic  noise  electronics.  Should 
this  prove  to  be  feasible,  as  is  expected  based  on  (lie  present  results,  it 
would  then  be  useful  to  consider  application  of  the  technique  up  to  the  100  krn 
region  either  by  use  of  a  more  intense  source  (several  kCi),  by  developing 
detection  techniques  allowing  a  lower  energy  threshold  (£  10-Z0  keV)  -  com¬ 
bined  with  a  low  energy  beta  source  like  Ni-oj,  oracoinbinationofthe.se  measures. 

Inasmuch  as  the  U.S.  Government  has  a  considerable  investment 
in  this  sonde  at  this  time,  we  believe  that  it  should  be  utilized  in  the  flight 
program  recommended  above. 
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